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Abstract. Cosmic rays extensive air sliowers (EAS) 
are associated with transient radio emission, wliicli 
could provide an efficient new detection method 
of high energy cosmic rays, combining a calori- 
metric measurement with a high duty cycle. The 
CODALEMA experiment, installed at the Radio 
Observatory in Nan^ay, France, is investigating this 
phenomenon in the 10^^ eV region. One challenging 
point is the understanding of the radio emission 
mechanism. A first observation indicating a linear 
relation between the electric field produced and the 
cross product of the shower axis with the geomag- 
netic field direction has been presented (B. Revenu, 
this conference). We will present here other strong 
evidences for this linear relationship, and some hints 
on its physical origin. 
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I. Introduction 

Radio detection is a promising detection technique 
of extensive air shower as it could provide at low- 
cost a calorimetric measurement of the shower with 
a high duty cycle. The first exploration phase of this 
technique is described in Allan's review of 1971 |]T|. 
The observations of the electric field in the East- West 
polarization are there summarized with the following 
equation of the electric field produced on the ground 
by an EAS: 
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Ep stands for the primary cosmic ray energy in eV, a is 
the angle between the shower axis and the geomagnetic 
axis, 6 the zenith angle and R radial distance to the 
shower axis. The value Rq is estimated for 6* < 35 and 
at two frequencies: 100±10 m at 55 MHz and ^140 m 
at 32 MHz. 

The radio detection technique has been neglected for 
a while in favor of the fluorescence and surface detectors 
developments. Nowadays, progress of fast and low cost 
electronics offers a new insight for radio detection. The 
experiments CODALEMA in France and LOPES in 
Germany are currently dedicated to the study of this 
technique in the lO^^ eV region [[21, 13, along with 
other preliminary experiments in Argentina H, fSl. This 
experimental work is supported by a renewal in theories, 
with several emission models and implementations such 



as the geosynchrotron emission, the boosted Coulomb 
field or the transverse current ll6l. llTll. ill. ll9l. ifTOl. 

II. The v X B polarization vector 

The synchrotron emission of an electron moving 
through the Earth's magnetic field can be calculated 
using the general formula of an accelerated relativistic 
particle, with v representing here the velocity of the 
particle and n the direction of the observer: 
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The second term of the equation E2 represents the 
geosynchrotron emission when v — — v x B is the 
Lorentz acceleration in the magnetic field B. In the 
special case of an observation point lying on the particle 
motion axis, we have v |[ n: 

E2 oc n X {(n — v) X (v X B)} 

oc {n.(v X B)} (n- v) - {n.(n- v)}(v X B) 

=0 on the axis 

cx -(1 - ■(;)(v X B) 

cx -v X B (3) 

thus the synchrotron electric field produced near the 
particle motion axis is at first order proportional to this 
cross product. 

The transposition can be made to an EAS, where v 
represents from now on the direction of the shower axis, 
and B corresponds to the geomagnetic field. The E oc 
—V X B relation is conserved in the case of a realistic 
EAS, as checked with a full Monte Carlo simulation of 
the geosynchrotron emission in a modified version of 
the program AIRES fTS), lfT4l . These simulation results 
are presented in another contribution to this conference 
(C. Riviere ifTSl ). Other emission models may give the 
same dependence. 

We can note that the modulus of the v x B vector 
is sin a, as in Eq. [T] However, the sin a corresponds 
here to the total field amplitude, not to the East- West 
polarization. The East- West component of the fields is 
then given by the projection of the v x B vector on the 
EW axis. 

III. The CODALEMA experiment 

The CODALEMA experimental setup, located in 
Nangay, France, is currently composed of an array of 
17 particle detectors overlapped with an array of 24 
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wide band active dipoles: 21 oriented along the East- 
West (EW) direction and 3 oriented along the North- 
South (NS) direction, disposed as shown on Fig. [T] The 
whole acquisition is triggered by the particle detector 
array, which also provides reference information on the 
shower characteristics. After offline filtering of the radio 
signals in the 23-83 MHz band, the radio information is 
reconstructed and then confronted to the surface detector 
data. Additional information about the CODALEMA 
setup and data analysis can be found in Ref. ifTTl . 



product V X B, we obtain the function plotted on Fig. [3] 
Making two simple assumptions, this can be expected 
to represent an event rate distribution: i) the electric 
field is proportional to the vector cross product, ii) the 
probability of detection is proportional to the electric 
field. This is realistic here because we are close to the 
detection threshold. 

We see that this simple model reproduces well the 
characteristics of the direction of the radio signals seen, 
in the two measured polarizations. 




Fig. \. Layout of the CODALEMA experimental setup. Scintillator 
stations are represented in pink squares, EW antennas in yellow, NS 
antennas in green. 

A first observation indicating a linear relation between 
the electric field produced and the cross product v x B 
of the shower axis v with the geomagnetic field di- 
rection B is presented in another contribution to this 
conference lfT2l . We will detail bellow evidences of this 
relation. The following analysis is based on data taken 
during 453 effective days of stable acquisition with this 
configuration of 24 antennas. A radio event is an event 
with at least three fired antennas in order to reconstruct 
a radio shower plane. A radio detected event is a radio 
event which matches the coincidence criteria with the 
ground detector (arrival direction difference < 20° and 
time difference < 200 ns). 

IV. Experimental sky coverage 

The arrival directions of the radio detected events 
present a strong anisotropy. The distributions of arrival 
directions are different when considering the EW or NS 
polarization, as shown on Fig. |2l EW oriented antennas 
tends to see more showers coming from the North, 
whereas the NS oriented antennas tends to see showers 
coming from East and West. 

The radio detection threshold of CODALEMA is 
around 10^^ eV at 50%. Let us consider the distribu- 
tion of events above this energjQ. If we multiply this 
distribution by the different projections of the cross 




Fig. 2. Arrival directions of the radio detected events containing at 
least one EW (left) or NS (right) polarized flagged signal. 
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Fig. 3. Product of CODALEMA particle detector acceptance above 
10^^ eV times the EW (left) or NS (right) component of the vector 
cross product (absolute values). 

Similar results are found in Argentina with a different 
geomagnetic field orientation, 35° to the North instead 
of 63° to the South in Nangay. On Fig. |4] left, are 
shown the arrival direction of the events detected with 
autonomous radio detectors [4J . These events have been 
self triggered with the radio signals of the EW polarized 
antennas, and are in coincidence with events seen with 
the surface detector of the Pierre Auger Observatory. On 
the right side of Fig. H] is plotted the absolute value of 
EW component of the vector cross product multiplied by 
an estimate of the Auger zenith acceptance distribution. 
The depletion of events on the North is also explained 
with the same model. 

V. Detection efficiency 

We consider here the events reconstructed with the 
array composed of the EW oriented antennas. The radio 
detection efficiency of the CODALEMA experiment is 
shown on Fig. |5] In gray, the detection efficiency is 
plotted against the energy estimated by the particle 
detector array, E. In black, the efficiency is plotted 
against E' = E.\{v x'B/vB)ew\, the energy multiplied 
by the EW component of the vector cross product, which 
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Fig. 4. Left: arrival directions of the self triggered radio events 
detected in Auger (left). Right: Auger surface detector acceptance 
estimation multiplied by l(v X S/vB)ew\ 



represents an estimation of the electric field amplitude. 
As expected, the efficiency rises faster in the second 
casell. 
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Fig. 5. Radio detection efficiency, as a function of the energy E 
(gray) or of the connected energy E' (black). 



Regardless the energy, the radio detection efficiency 
also increases with |(v x 'B/vB)ew\^ as shown on 
Fig. |6] for the events above 10^^ eV. The quasi linear 
tendency is in favor of the assumption made earlier that 
the detection efficiency is proportional to the electric 
field close to the energy thresholcj^. 

VI. Signal sign 

The sign of the components cross product — v x B 
varies with the EAS arrival direction. The dependences 
of the EW and NS components are shown on Fig.]?] The 
EW (resp. NS) component is positive on North (East) 
and negative on South (West). 

Experimentally, one antenna signal is defined as the 
extremum of a time signal filtered in the 23-83 MHz 
bancfl. The sign of the signal is simply the sign of 

-Width of a Fermi-Dirac fit of the efficiency curve: 0.171 (E) and 
0.137 (E'). Rq: 0.171 with a sina correction, and 0.157 with a (1 — 
cos a) correction. 

^The correlation is much weaker with sin a or (1 — cos a). 

'*Filtering of simulated signals indicates that the sign is conserved 
during the process, but anyway the sign flip with the arrival direction 
is more important than the absolute sign. 
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Fig. 6. Radio detection efficiency vs. the EW component of the vector 
cross product [(v X S/vB)e^y\, for the events above 10^'' eV. 




Fig. 7. Sign of the EW (left) and NS (right) components of the vector 
cross product — v X B. The scale is does not have importance here. 



this extremum. For an event with several antennas fired 
for the considered polarization, the sign of the event is 
chosen to the majority sign among the different signals. 
This reduces the influence of noise, and is of course 
more efficient in EW polarization than in NS as more 
antennas are available. The EW component of — v x B 
is also generally bigger than its NS component, thus the 
measurement of the sign is generally easier in the EW 
polarization. 

The experimental signs of the radio detected events 
are shown on Fig. [8] In the EW polarization, events 
coming from North generally have a positive sign and 
events coming have a negative sign. The apparently 
uniform distribution of the negative event is simply a 
statistic effect where appear the negative events coming 
from South plus a small fraction of the very numerous 
events coming from North. In the NS polarization, events 
coming from East are generally positive and events 
coming from West are generally negative. The overall 
distributions of the sign with the arrival direction is the 
same as the sign distribution of — v x B: the agreement 
is 93% for the EW polarization (with 19 antennas 
maximum), and 78% for the NS polarization (with 3 
antennas maximum). 

VII. Conclusion 

Important progress was achieved recently in the field 
of radio detection of high energy cosmic rays. A 
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Fig. 8. Experimental signs of the radio signals seen in Nanjay, with the EW (left) or NS (light) oriented antennas. Red dots represent positive 
signals and blue circles negative ones. 



complete understanding of the electric field production [15] C. Riviere, this conference (3P* ICRC, Lodz 2009) 

mechanism is necessary to go back from the radio sig- |a rXiv:0906.26981 . 

nals to the primary cosmic ray characteristics (direction, 

energy, nature). Different theoretic approaches of the 

radio emission are under investigation, some of which 

predict at first order a linear dependence of the radio 

electric field with the vector cross product — v x B. 

This dependence offers a simple interpretation of the 
experimental observations of the CODALEMA exper- 
iment in two polarizations, such as the anisotropy of 
detection efficiency close to the threshold or the signs of 
the signals in each polarization. This polarization factor 
is a major effect to take into account when observing 
cosmic rays with radio antennas. 

References 

[1] H.R. Allan, Radio emission from extensive air showers. Progress 
in Elementary Particles and Cosmic Ray Physics 10 (1971) 171. 
[2] D. Ardouin et al., Nucl. Instrum. Meth^ A 555 (2005) 148. 
[3] H. Falcke et al.. Nature 435 (2005) 313. 

[4] B. Revenu for the Auger Collaboration, Third International 

Workshop on the Acoustic and Radio EeV Neutrina detection 

Activities ARENA, Rome, 2008, NIMA (in press). 
[5] J. Coppens for the Auger Collaboration, Third International 

Workshop on the Acoustic and Radio EeV Neutrina detection 

Activities ARENA, Rome, 2008, NIMA (in press). 
[6] H. Falcke and R Gorham, Astroparticle Physics 19 (2003) 477. 
[7] T. Huege et al., Astroparticle Physics 27 (2007) 392. 
[8] T. Gousset, J. Lamblin and S. Valcares, Astroparticle Physics 31 

(2009) 52. 

[9] N. Meyer- Vemet, A. Lecacheux and D. Ardouin, Astronomy and 

Astrophysics 480 (2008) 15. 
[10] K. Werner and O. Scholten, Astroparticle Physics 29 (2008) 393. 
[11] D. Ardouin and the CODALEMA Collaboration, Astroparticle 

Physics (2009). 

[12] B. Revenu for the CODALEMA Collaboration, this conference 

(31"* ICRC, Lodz 2009) | arXiv:0906.2832 1. 
[13] J.S. Sciutto, AIRES, http://www.fisica.unlp.edu.ar/auger/aires/ 
[14] M.A. DuVernois, B. Cai and D. Kleckner, Proc. of the 29th 

ICRC, Pune, India (2005). 



